Application No. 10/534,735 PATENT 
Response Dated 03/30/2008 

Reply to Office Action of 12/30/2008 Agent's Docket No. 2173-230 

Amendments to the Specification: 

Please replace paragraph [0050] of the published application with the following corresponding 
amended paragraph: 

[0050] Matlab Simulation Listing 

% viterbi demodulator for phase excursion 

% 

% number of state transitions 
K = 2000; 

% number of states 
M = 30; 

% phase angle of each state 

state_angle = ( 2* . quadrature . pi/M) * [0:M-1]'; 

% model the actual signal phase 
phi_sig = .03 * [1:K] ' ; 

% actual signal 
SNRdB = -6; 
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SNR = 10 A { 0 . l*SNRdB) ; 
u = exp (lj*phi_sig) ; 

x = randn(K, 1) + 1 j *randn (K, 1) + sqrt(SNR) * u; 
phi_meas = angle (x); 

% metric for phase state 

% at each update there is only one surviving path for each state 
Mt = zeros (M, 1 ) ; 

% paths, and initiate the k=l state; equal probability that we 
can start from 

% any state 

P = zeros (M, K) ; 

P(:,l) = [1:M] '; 

Pnew - P; 

% loop over the overall set of K state transitions 
I = ones (M, 1) ; 

for k=2:K 

% determine the path metric of the new segment from state 
k-1 to state k 

gl = phi_meas(k) * ones{M,l) - state_angle; 
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g2 = gl + 2* ■ quadrature .p i*I; 
g3 = gl - 2 * -^- quadrature . pi*I ; 
M_new = (min(abs ( [gl,g2,g3] ') ) ' ) - A 2; 

% determine the new path for state m 
for m=l :M 

if m == 1 

mml = M; 
mO = 1; 
mpl = 2; 
elseif m == M 

mml = M-l ; 
mO = M; 
mpl = 1; 
else mml = m-l; 
mO = m; 
mpl = m+1; 

end 



[gl,i] = min ( [Mt (mml) ,Mt (mO) ,Mt (mpl) ] ) ; 
if i == 1 

Pnew(m, l:k) = [P {mml, 1 : k-1) ,m] ; 

Mt (m) = Mt(mml} + M_new (m) ; 
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elseif i ==2 

Pnew (m, 1: k) = [P {m0, 1 : k-1) ,m] ; 
Mt(m) = Mt(m0) + M_new(m) ; 

else 

Pnew (m, 1 : k) = [P (mpl, 1 : k-1) # m] ; 
Mt (m) = Mt(mpl) + M_new(m); 

end 

end 
P = Pnew; 

end % end of the K transition periods 

figure (1) ; plot (P' ) ; 
hold on; 

[gl, i] = min (Mt ) ; 
plot (P (i, :}, 'or'); 
hold off; 

% unwrap P 

Puw = zeros (K, 1 ) ; 

Puw(l) = P(i,l) ; 

offset = 0; 

for k=2:K 
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Puw(k) = P(i,k) + offset; 
if Puw(k) - Puw(k-l) > 1 

offset = offset - M; 

Puw(k) = Puw(k) - M; 

end 

if Puw(k) - Puw(k-l) < -1 
offset = offset + M; 
Puw(k) = Puw(k) + M; 

end 

if abs(Puw(k) - Puw(k-l)) > 1 

error {' discontinuity in Puw'); 

end 

end 

figure (2) ; 
plot (Puw) ; 

figure (3) ; 
plot(P(i, :) , 'r') ; 
hold on; 

plot (mod (phi_sig * M/(2 * . quadrature . pi ) , M) , ' g ' ) ; 
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hold off; 
grid on; 

% calculate the integration loss 
Loss = 20 * logl0(abs{u' * exp (1 j *state_angle (P (i, : ) ) ) ) /K) 
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